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Abstract. Symmetry arguments are used to develop a spin Hamiltonian for the 
description of the complex magnetic ordering in HoMnOa. Using a novel application of 
the Landau Lifshitz Gilbert dynamic torque equations to this model of the frustrated 
Mn ions on an AB stacked triangular antiferromagnetic, it is shown that the four 
principal spin configurations observed in this compound are stabilized. Ho-Mn coupling 
is found to be a consequence of an unusual trigonal anisotropy term which is responsible 
for simultaneous Mn spin reorientation and onset of Ho magnetic order. Based on 
these microscopic considerations, a mean-field Landau-type free energy is derived which 
reproduces the succession of observed temperature driven magnetic phase transitions 
at zero field, including re-entrant behavior. In addition, our analysis suggests that the 
basal-plane magnetic order should be slighty incommensurate with the lattice. 
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Hexagonal HoMnOs is one of the most studied of the rare-earth manganites and has 
a magnetic field - temperature phase diagram that exhibits a multitude of complex spin 
structures [H El El lU El E] . Key features of the magnetic phases in this stacked triangular 
antiferromagnet include a series of reorientation transitions (with re-entrant behavior) 
at about 5 K and 40 K involving Ai^-stacked triangular layers of Mn spins which form 
the familiar basal-plane 120° spin structure below the Neel temperature Tjy = 76 K. 
Low temperature c-axis spin alignment of AA-stacked Ho ions is also observed which is 
concomitant with Mn spin reorientation. The nature of the coupling between Mn and 
Ho ordering has remained puzzling for close to two decades. 

In this Communication, a formulation of the spin Hamiltonian appropriate for a 
description of the magnetic phases of HoMnOs is developed on the basis of symmetry 
arguments. The principal Mn spin configurations observed in the magnetic phase 
diagram are found to depend on the signs of the coefficients of in-plane sixth-order 
anisotropy and inter-plane exchange coupling. These states are determined through 
a numerical solution of the dynamic Landau-Lifshitz-Gilbert (LLG) equation which 
includes finite temperature effects. Theoretical arguments proposed over twenty years 
ago that even weak inter-layer coupling of AB stacked triangular antiferromagnetic 
layers leads to an incommensurate distortion of the usual 120'' period-3 spin structure 
[7] are confirmed by these simulations. The crucial role of Ho-ion ordering in driving the 
transitions between the various Mn-ion ordered states is also explicitly demonstrated 
|8]. Since strong planar anisotropy in the case of Mn spins (S _L c) and axial anisotropy 
in the case of Ho spins (So || c || z) leads to S ■ So = 0, the usual exchange coupling 
mechanism should be absent in equilibrium (but can contribute to magnetic excitations 
[9]). However, the hexagonal crystal symmetry of this compound contains a 6- fold screw 
axis which allows for a fourth-order trigonal anisotropy term of the form [10] 



where the sum is over near neighbors. Such a term can occur provided that either 
Ho or Mn spin configurations are antiferromagnetic between AB layers, but not both. 
The consequences of this interaction are examined within the framework of a simple 
Landau-type free energy derived within a mean-field approximation. This demonstrates 
that the trigonal coupling term causes the series of Mn-spin reorientation transitions 
observed in zero applied magnetic field. We note that a phenomenological Landau-type 
free energy based on group theoretic arguments was previously developed to explore the 
magnetic phase diagrams of hexagonal RMnOs compounds, which is complementary to 
the present work [T2] . 

Our investigation may be compared with other proposals regarding microscopic 
mechanisms of coupling between rare-earth and Mn ions in the hexagonal manganites. 
A recent discussion of Mn-Yb coupling in the sister compound YbMnOs (which does not 
exhibit zero-field Mn-reorientation transitions) suggests that it could be due to dipolar 
or Dzyaloshinksii-Moriya (DM) interactions [13]. Dipolar effects were estimated to be 
an order of magnitude too small and either effect would require that the rare-earth ions 
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order at the same temperature as the Mn ions (since both types of interactions involve 
a hnear couphng between the two types of spins), which is not the case in HoMnOs. 

An intriguing alternative explanation has been put forth which suggests that the 
slight distortion of the perfect triangular symmetry can give rise to asymmetric inter- 
layer exchange interactions [I2l [H] . A change in sign of this exchange coupling occurs 
in response to the temperature-driven lattice distortion which in turn drives the Mn- 
spin reorientations. The significance of such distortion-driven interlayer-coupling may 
be considered secondary in view of the observation made here that ordinary inter- 
layer exchange should indeed be present (since the in-plane Mn modulation should 
not be exactly period-3). In addition, it is well established that there are simultaneous 
changes in the Ho magnetic order at the Mn-spin reorientation transitions which must 
be accounted for in a complete microscopic model. We do show below, however, that 
the distortion of the perfect triangular lattice resulting in anisotropic in-plane exchange 
interactions serves to enhance the incommensurability of the Mn spin modulation (in 
addition to the ordinary exchange coupling between AB stacked layers mentioned above) 

Terms which appear in the spin Hamiltonian (or Landau-type free energy) of a 
magnetic system can be constructed with the requirement of invariance with respect 
to the generators of the crystal space group in the paramagnetic regime. In the case 
of hexagonal RMnOs compounds, the crystal symmetry group is non-symorphic PG^cm 
with generators given by a screw rotation {Cg'lOOl} and glide plane {o"^|00|} |TT]. Terms 
involving lattice vectors (such as dipole contributions) are omitted here for simplicity. 
In the present case, allowed anisotropic terms at second order are the axial/planar type, 
{S^y, as well as the DM interaction. The lowest order term which gives rise to a 
dependence of the orientation of planar spins relative to the crystal axes occurs at sixth 
order. 

In the remainder of this Communication results are presented on LLG simulations 
involving only Mn spins and an analysis of a mean-field Landau-type free energy with 
Ho-Mn trigonal coupling. We conclude with a discussion of incommensurabiltiy arising 
from distortions of perfect triangle. 

With the goal exploring a simplified spin Hamiltonian involving only the Mn ions 
we consider the following contributions 

+ ^^[(•S'l + ^S'lf + (SJ - iS'lf]. (2) 

3 

The first term includes only isotropic in-plane near-neighbor exchange interaction of the 
perfect triangular lattice, J = 1 as well as near-neighbor exchange between A and B 
triangular layers, J' . Interlayer exchange is assumed here to be weak (as in YMnOa 
[I6]). As noted previously [7], this contribution is zero in the case of the 120° spin 
structure since the coupling between a spin on layer A and its three neighbors on layer 
B appears as J'S^ • (Sf -f Sf + Sf ). Strong planar anisotropy is assumed to arise from 
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a large single- ion term with D < 0. The sign of the in-plane anisotropy coefficient E 
determines the orientation of S relative to the crystal a-axis in a temperature regime 
below Tat. 

This model Hamiltonian was used to demonstrate the four principal Mn spin 
configurations identified in RMnOs compounds. These are conveniently illustrated in 
Fig. 7. of Ref. [13] and are labeled ri^4 and can be characterized by two features: (1) 
inter-plane ferromagnetic or antiferromagnetic spin configurations, and (2) intra-plane 
spin orientations being parallel or perpendicular to a basal plane a-axis. 

Spin configurations were determined numerically as a solution to the dynamic LLG 
equations, which may be expressed as the the familiar torque equation, involving the 
gyromagnetic ratio 7, plus a damping term involving the parameter a, as 

f = -rT^SxH.„-^Sx(SxH.„). (3) 

The LLG equations have proven very useful in the determination of subtle differences 
in equilibrium magnetic structures which involve long-range magnetostatic effects in 
thin films but have occasionally also been applied to spin systems with only short-range 
exchange interactions pTl [18] . They are well suited for the study of highly frustrated 
magnetic systems and may be viewed as an alternative to Monte Carlo simulations. 
The effective field is given by He// = —51-L/5S + Yith where Yith is the stochastic 
field describing thermal fluctuations within the framework of Langevin dynamics. Spin 
vectors were located on interpenetrating 12 x 12 x 12 AB stacked triangular lattices. An 
Euler numerical integration scheme with adaptive time steps was employed using 7 = 1, 
a = 0.1 where equilibrium spin structures are determined by averaging the long-time 
behavior of Si(t). 

The Neel order parameter can be calculated as the root mean square of the three 
sublattice magnetizations as in Ref. [18]. For simplicity, a large single ion anisotropy 
D = — 1.0 was assigned (although larger than the experimentally determined value [3], 
this difference is not relevant to the present calculation which simply demonstrates the 
four zero-field Mn spin states) and small interlayer exchange \J'\ = 0.01 and small 
in-plane anisotropy \E\ = 0.01, we find T/v — 0.43. 

The four principal Mn spin phases result from these simulations by changing the 
signs of J' and E. An example spin configuration determined at T = 0.40 is shown in 
Fig. 1 which correspond to the Fi phase (S^ = —Sb, S ± a). Note that the effect 
of thermal fluctuations is to induce deviations from the idealized structure. The other 
three principal phases were also observed in the numerical results, dependent on the 
signs of J' and E , as follows: F2 with J' > 0, E < 0; F3 with J' < 0, E < 0; F4 with 
J' < 0, ^ > 0. 

We note that the near-neighbor exchange interaction between Mn spins in HoMnOa 
has been estimated to be J = 2.44 meV ~ 28 from spin- wave dispersion curves [3]. 
Using s = 2 associated with Mn^"*", we find T/v — 0.43Js^ ~ 112K, considerably larger 
than the observed 76 K. The reasons for this discrepancy are not clear but may be due 
to the assumed value of J'. 
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Figure 1. LLG simulation results using J = 1.0, J' = 0.01, D = -1.0, = 0.01 
T ~ 0.4. This corresponds to the case of phase Fi. An outline of the unit cell is shown. 



The crucial role of the Ho-ion spins in driving Mn-spin reorientation transitions 
is made clear by considering a simple Landau-type free energy, derived from mean 
field theory, which accounts for the microscopic interactions identified above. This 
approach has the advantage of revealing explicitly the competition between the different 
microscopic energy and entropy [S) contributions through the relation F = E — TS, 
where thermal effects are calculated in a manner similar to the Brillioun function [191 120] . 
Single-ion anisotropy associated with the Ho ions is assumed to be positive {Dq > 0) 
giving So II c || z. Keeping only relevant terms to sixth order yields (in units of the 
Boltzmann constant): 



where A = a{T — Tjy) and Aq = a^iT — T^q) B = bT, Bq = b^T, C = cT, Cq = CqT. 
Here, the Mn^+ spins order at the Neel temperature T/v = — J(Q)s^/a. For the Ho^"*" 
ions, the corresponding parameter is T^o = (— ^o(Qo) + Dq)]"^ /oq which involves the 
total angular moment [1] j = 8. Expressions for the entropy constants are given in 
Ref.[20] and yield a = 2, 6 = 0.867, c = 0.610, ao = 2.667, 6o = 1-43, and cq = 1.21. 
The remaining model parameters involving the exchange and anisotropy are estimated 
to reproduce the approximate sequence of phase transitions observed at zero field. We 
note that this approach is expected to be increasing less quantitative as the temperature 
is lowered where the Landau assumption of a small order parameter S breaks down. 

For spin structures with a periodicity of 3 relative to near-neighbor spacing on a 
triangular lattice, J(Q) = — 3Ji (also see below). For this analysis, we set = 76K 
which yields the mean-field estimate Ji = 1.15meV, more than a factor of two smaller 
than the measured value. This difference is consistent with the expected role of spin 
fluctuations in quasi-2D frustrated antiferromagnets (a suppression of T/v). An estimate 
of T/vo can be made using the measured Curie- Weiss behavior fT], which yields the 
estimate Tq = —UK. This value may be compared with the above expression for 
Tatq using (5 = and Jo(0) = 6Joi where Jio is the near-neighbor exchange interaction 
between Ho ions in the triangular plane. For simplicity, we ignore the anisotropy Dq so 
that the estimate Jio — O.li^T can be made. This result is consistent with very weak Ho- 




(4) 



A model of magnetic order in hexagonal HoMnO^ 



6 



Ho spin coupling. Assuming a period-3 near-neighbor modulation of the Ho moments 
then gives T/vo — 87^, associated with a pure Ho system. 

In the free energy (4), represents the angle between Mn spins and the triangular 
lattice defined by the Ho ions. The trigonal i^-term is minimized by = nvr/S, 
independent of the sign of whereas the basal-plane anisotropy ii^-term favors 
= n7r/3 for E < and (2n + l)7r/6 for E > 0. 

Numerical minimization of F{S, Sq, (p) reveals a variety of phase transition 
sequences as the temperature is lowered below the paramagnetic phase {S = Sq = 0), 
dependent on the assumed values for E and K. Variations occur mainly through the 
sign of E and the magnitude of K. For example, the using K = 25 K and E = 1 K (for 
comparison note that 6T/v — 66 K) yields the results shown in Fig. 2. The first ordered 
state is characterized hy S 0, So = and = (2n + l)vr/6, corresponding to r4 or 
Fi. In the narrow range between around T ~ 37 we find that decreases rapidly to 
zero (0 = nvr/S) with the concomitant development of Ho-spin order. Analysis of the 
free energy reveals that the Ho ordering occurs at Tho = Tjsfo + K'^/ {SEoq). The narrow 
region around 37 K (where < < 30°) corresponds to phase F5 or Fg. In the range 
T 37 K to T c:^ 12 K, the order corresponds to F2 or F3. Finally, below T ^ 12 ii", 
becomes non-zero again, corresponding to re-entrance to F5 or Fg. 

The trigonal term is allowed by symmetry only if the magnetic order alternates in 
sign along the c-axis for either Mn or Ho spins (but not both) and is ferromagnetic for 
the other species. A conclusion of Ref. ^ is that Ho order is antiferromagnetic between 
triangular planes so that for phases where 5*0 7^ in the present analysis, the Mn order 
must be ferromagnetic along the c-axis. This suggests that the main sequence of phases 
in HoMnOs as the temperature increases is given by Fq — ?■ F3 — )■ F4, with a narrow region 
of stability of the intermediate phase Fg between F3 and F4. This complex sequence of 
transitions, involving re-entrance of the state Fg, is consistent with experimental results 
(see, e.g., Ref. 

Other choices of parameter values yields different phase transition scenarios. For 
example, a smaller value of K = 10 K stabilizes the sequence Fg — > F4, with increasing 
temperature. With E < 0, (j) = m:/3 is stabilized (F2 or F3), which yields simultaneous 
Mn and R ordering at Tjy. These scenarios have been reported to occur in ErMnOs, 
YbMnOg and TmMnOg. 

The pricipal impact of a magnetic field applied along the c axis is the suppression 
of the F3 phase (0 = r;,7r/3). A simple analysis illustrates how this can arise within the 
Landau free energy model. An applied field will induce a uniform contribution to the 
spin density giving 

s(r) = m + [Se'^" + So*e-'^''" + c.c] (5) 

and the free energy functional F[s(r)] (with a Zeeman term added) expanded in powers 
of s will contain numerous additional terms involving the magnetization m = xHz 
Some of the additional terms serve to effectively introduce a field dependence to the 
coefficients which appear in the free energy (4). For example, both anisotropy terms are 
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Figure 2. Results of numerical minimization of the mean- field free energy (4). 

renormalized as — > Kh = K + K^m? and E — )■ Eh = E + Eqiti? where Kq and Eq 
are new parameters. The stabihty of the phase from the numerical minimization of 
the free energy (4) is found to be very sensitive to the values of both Kh and Eh- For 
example, with the nominal zero-field value Eh = 1-0, this phase is stable in the region 
12 K < T < 37 K. A field induced increase to a value Eh = 1.1 would shrink this to 
a temperature interval 16 K < T < 34 K and at Eh = 1.3 the r3 phase is suppressed 
entirely in favor of the r4 phase. Similar sensitivity is found with variations in Kh- 
We note that a term of the form mS"^ cos(30) can also occur but only if the Mn spins 
order anti-ferromagnetically between planes. This term thus enhances the stability of 
the Fi phase. These results are consistent with experimental magnetic phase diagrams 
of HoMnOs where above a relatively small applied field of 3-4 T the F4 spin structure 
occupies most of the phase diagram and the F3 state is absent. In addition, the Fi phase 
is observed at low temperatures but only at fields above 1.5 T. At lower temperatures, 
where the Landau model assumptions are less reliable, we find that the intermediate 
phase Fg is stabilized. A wide variety of scenerious have been proposed on the nature 
of the spin ordering in this region of the phase diagram based on experimental work. 

The above analysis was made with the assumption that the spin orderings are 
commensurate with the lattice. Due to both small inter-plane coupling of Mn spins as 
well as the small distortion of the triangular lattice, we argue that the magnetic ordering 
should in fact be slightly incommensurate. We confirmed in the numerical solutions to 
the LLG equations the expected incommensurabilty of the in-plane modulation driven 
by interlayer coupling J' of AB stacked triangular layers [7]. This was observed as 
a deviation from 120*^ in the average Mn inter-spin angle as a function of small J'. 
The effect of the distortion of the triangular lattice (characterized by the deviation 
of the parameter x from 1/3) is seen by considering the wave vector which minimizes 
the Fourier transform of the in-plane Mn exchange interactions, J(Q). Referring to 
Ref. [1], vectors connecting site to its nearest neighbors are given by rio = (|a;, x), 
r20 = (fa;, -x), rgo = (fx - i,x - 1), = (fx - 1,-x), rgo = (fx - l,x), 
Tqq = (|x — |, 1 — x) relative to crystallographic vectors (ax, 6y) where b = \/3/2. In 
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this description, the simple hexagonal lattice is defined by the Ho-ion sites. There are 
two Mn-Mn bond lengths, {\/3x^)a and (a/Sx^ — 3a; + l)a, corresponding to exchange 
interactions [16], Ji > J2, giving 

3 

J(q) = 2 Ji cos{-xqx) cos{qyx) 

+ 2J2[cos(^xg^ - ^g^.) cos{qyX - Qy) 
3 

+ cos(-xg^ - q^) cos(gyx)] (6) 

where {qx,qy) = {(iQx,bQy)- With J2 = Ji and x = |, this function is minimized 
by near-neighbor spins having a relative angle of 120*^. Deviations from this value as 
a function of either x (with J2 = Ji) or J2/J1 (with a;=0.322, the observed value in 
HoMnOs pj) are shown in Fig. 3. 




Figure 3. Results of numerical minimization of J(Q) (4) showing the impact of the 
angle between near- neighbor Mn spins due to the lattice distortion x (with J2/J1 = 1) 
and also due to J2/J1 7^ 1 (at the value 2:=0.322). 

In addition to demonstrating that the four principal Mn-spin states in HoMnOa 
are a consequence of a simple spin Hamiltonian, our mean-field analysis has also shown 
that an unusual trigonal anisotropy term leads to Mn-Ho couplng which accounts for the 
observed complex phase transition sequence at zero field. This investigation has revealed 
explicitly the relation between Ho spin order and Mn spin reorientation. Extensions of 
the LLG approach, which is not limited by the same assumptions as Landau theory, to 
include this coupling and an applied magnetic field to examine phase diagrams and spin 
excitations are planned. Our work has also provided strong theoretical arguments which 
support a scenario where the in-plane modulation characterizing Mn spin order in many 
RMnOa compounds should exhibit a slight incommensurability away from the normally 
assumed period-3 modulation. This deviation from a perfect 120° spin structure has 
not yet been reported, likely a consequence of the smallness of the incommensurability, 
but deserves a more detailed investigation. 

We thank J. Mercer, G. Quirion and S. Curnoe for enlightening discussions. This 
work was supported by NSERC of Canada and the Atlantic Computational Excellence 
Network (ACEnet). 



A model of magnetic order in hexagonal HoMnO^ 



9 



References 

[I] Munoz A, Alonso J A, Martinez-Lope M J, Casdis M T, Martinex J L and Ferndndez-Diaz M T 

2001 Chem. Mater. 13 1497 
[2] Fiebig M, Dcgcnhardt C and Pisarcv R V 2002 J. Appl. Phys. 91 8867 

[3] Vajk O P, Kenzelmann M, Lynn J W, Kim S B and Cheong S-W 2005 Phys. Rev. Lett. 94 087601 
[4] Lorenz B, Yen F, Gospodinov M M and Chu C W 2005 Phys. Rev. B 71 014438 
[5] Yen F, dcla Cruz C, Lorenz B, Galstyan E, Sun Y Y, Gospodinov M and Chu C W 2007 J. Mater. 
Res. 22 2163 

[6] Camirand J, Poirier M, Loreynne P.-G. and Revcolevschi A 2009 Phys. Rev. B 79 094423 
[7] Sakakibara T 1984 J. Phys. Soc. Japan 53 3607 

[8] Nandi S, Krcyssig A, Tan L, Kim J W, Yan J Q, Lang J C, Haskel D, McQueeney R J, and Goldman 

A 1 2008 Phys. Rev. Lett. 100 217201 
[9] Talbayev D, LaForge A D, Trugman S A, Hur N, Taylor A J, Averitt R D, and Basov D N 2008 

Phys. Rev. Lett. 101 247601 
[10] Simpson J A and McMorrow D F 1995 Phys. Rev. B 51 16073 

[II] Plumer M L 2008 Phys. Rev. B 76 144411; ihid 2008 78 094402 

[12] Munawar I and Curnoe S H 2006 J. Phys.: Condens. Matter 18 9575 

[13] Fabreges X, Mirebeau I, Bonville P, Petit S, Lebras-Jasmin G, Forget A, Andre G and Pailhes S, 

2008 Phys. Rev. B 78 214422 
[14] Fabreges X, Petit S, Mirebeau I, Pailhes S, Pinsard L, Forget A, Fernandez-Diaz M T and Porcher 

F 2009 Phys. Rev. Lett. 103 067204 
[15] Zhitomirsky M E 1996 Phys. Rev. B 54 353 

[16] Sato T J, Lee S-H, Katsufuji T, Masaki M, Park S, Copley J R D and Takagi H 2003 Phys. Rev. 
B 68 014432 

[17] Feng X and Visschcr P B 2001 J. Appl. Phys. 89 6988 

[18] Szunyogh L, Lazarovits B, Udvardi L, Jackson J, and Nowak U 2009 Phys. Rev. B 79 020403(R) 
[19] Bak P and von Boehm J 1980 Phys. Rev. B 21 5297 
[20] Plumer M L 1991 Phys. Rev. B 44 12376 



